Strongly coupled quantum criticality with a Fermi surface in two dimensions: 
fractionalization of spin and charge collective modes 
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We describe two dimensional models with a metallic Fermi surface which display quantum phase 
transitions controlled by strongly interacting critical field theories below their upper critical di- 
mension. The primary examples involve transitions with a topological order parameter associated 
with dislocations in collinear spin density wave ("stripe") correlations: the suppression of disloca- 
tions leads to a fractionalization of spin and charge collective modes, and this transition has been 
proposed as a candidate for the cuprates near optimal doping. The coupling between the order 
parameter and long-wavelength volume and shape deformations of the Fermi surface is analyzed 
by the renormalization group, and a runaway flow to a non-perturbative regime is found in most 
cases. A phenomenological scaling analysis of simple observable properties of possible second order 
quantum critical points is presented, with results quite similar to those near quantum spin glass 
transitions and to phenomenological forms proposed by Schroder et al. (Nature 407, 351 (2000)). 
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I. INTRODUCTION 

An important property of most quantum phase tran- 
sitions studied to date in systems with a metallic Fermi 
surface in spatial dimensions d > 2 is that the critical 
field theory for the order parameter is a free Gaussian 
theoryoa. This result has its origin in the fact that 
the order parameters considered can be expressed as a 
fcrmion bilincars, and consequently the order parameter 
fluctuations are efficiently overdamped and suppressed 
by fermionic particle and hole excitations near the Fermi 
surface. The temperature dependencies of physical ob- 
servables near the qua^um.£^ical_point have been per- 
turbatively computedBM OiMEJ, .-and these can be 
understood as. corrections to scaling^ at the Gaussian 
critical point£3. 

The quantum and thermal fluctuations near quantum 
phase transitions in Fermi systems have been used to 
interpret experiments in a variety of correlated electron 
materials. In the cuprate superconductors, the anoma- 
lous properties in the normal state have been described 
in terms of the proximity to a quantum phase tran- 
sition associated with the onset of spin density wave 
(SDWXpcpGiiarge density wave (CDW) order in a Fermi 
liquidE^ota. However, the observed anomalous proper- 
ties extend to quite high temperatures, and it would be 
preferable to explain these in theories with stronger non- 
linearities among the order parameter modes. Further, 
many of the anomalous properties extend to the optimal 
doping regime where there is no strong evidence of a large 
correlation length of the SDW and CDW orders. 

This paper describes an alternative class of quantum 
critical points whose critical theories are non-Gaussian, 
and which remain strongly coupled even in the presence 
of a Fermi surface because they are below their upper 
critical dimension. Only in such theories can the order 
parameter relaxation rate, and possibly the quasiparticle 
energy width, be genetically equal to ksT times a univer- 
sal numerical constants (T is the absolute temperature). 



Examples of such quantum phasg transitions abound in 
insulators and superconductorsEZBlij, and some of these 
have been used to explain loAfcJmanecature properties of 
the cuprate superconductorsl3E^K3Ei This paper will 
discuss transition in metals with a Fermi surface, and so 
the possible regime of applicability is restricted to higher 
temperatures where superconductivity is absent. 

Our primary focus will be on a transition with a non- 
local "topological" order associated with certain defects 
in the SDW/CDW order. However, long-range SDW or 
CDW order will not be present on either side of the criti- 
cal point. This transition, (shown in Fig. ^| below) is being 
offered as a candidateEjlHi-|fpij a possible optimal dop- 
ing quantum critical pointE§Eil in the cuprate supercon- 
ductors. This proposal for the optimal doping quantum 
critical point must be distinguished from a conventional 
SDW ordering transit, ion in a superconductor discussed 
recently by one of usOEj; the latter transition occurs at 
lower doping concentrations (near 1/8) and was used to 
predict and explain a number of recent low temperature 
neutron scattering and STM experiments. Here we are 
interested in higher dopings and temperatures, where, as 
we have already noted, the correlation length of conven- 
tional SDW/CDW order is surely quite small, but, as we 
shall discuss below, correlations in a topological order 
may be far more robust. 

We also hope for applications of our theory in other 
correlated materials, like the heavy fermion systems, 
where many experiments are in dis agree ment with the 
conventional weakly coupled theories 2 IrT In particular, 
we will show that the strongly-coupled critical points pre- 
dict a scaling structure for physical observables whiph is 
quite similar to those proposed by Schroder et alEB Re- 
lated scaling structures are also known to appear near 
strongly-coupled spin glass quantum critical points, and 
this is reviewed in Appendix |a|. 

There is one simple route to obtaining a strongly cou- 
pled field theory in the presence of a Fermi surface that 
has been discussed earlierEI: use restrictions from momen- 
tum conservation to prevent the order parameter from 
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coupling linearly to low energy fermionic excitations. For 
example, a transition with the onset of SDW order at 
the wavevector K, but the geometry of the Fermi sur- 
face such that no two points on the Fermi surface are 
separated by K. In such a situation the order parameter 
excitations remain undamped, and the fermions appear 
to be innocuous bystanders to the transition. However, it 
turns out that the Fermi surface can still be relevant for 
the structure of the critical theory in certain situations; 
these effects are related to those that arise in the transi- 
tions we do discuss below, and so we defer discussion of 
this point until Section |y|. 

Turning to transitions characterized by a non-local, 
'topological' order parameter, a situation that has been 
much discussed in the recent literatureE3 has the electron, 
cL (5 =f)J. is a spin index) fractionalize into a spinlcss 
charge e boson, b, and a neutral spin 1/2 fermion f q : 
c\ = 6/J This fractionalization transition is described by 
a Zi gauge theoryEH, and both b and / ; carry a unit gauge 
charge; the physical fermion, cj, is, of course, gauge neu- 
tral. If we denote the Zi gauge field by cr^ (i, j are site la- 
bels), then terms like aijtrtfji will generically appear in 
the effective HamiltonianEI In the fractionalized phase, 
where the fluctuations of are effectively quenched, 
this term is a BCS-like pairing of the / s fermions, and it 
implies that the Fermi surface will generically be gapped 
out (except possibly at special points if the pairing is 
anisotropic). Consequently, this fractionalization transi- 
tion does not provide us with a candidate quantum phase 
transition in the presence of a Fermi surface, and we will 
not consider it further in this paper. 

The main focus of this paper shall be on a fractional- 
ization transition that does less damage to the integrity 
of the electron: the electron (and hence its Fermi surface) 
remains intact on both sides of the transition, but charge 
neutral collective modes (in the particle-hole sector) do 
fractionalize; only at the quantum critical point, and in 
the associated non-zero temperature quantum critical re- 
gion, is the electronic quasiparticle ill-defined. Zaanen 
et al.c5 have recently given an appealing pictorial de- 
scription of such a transition. A transition with closely 
related physical content was introduced in Ref. |23|; we 
will follow the latter approach here, and will describe the 
physical content and experimental motivation below in 
Section ||. There have also been earlier discussions of 
the fractionalization of order parameters in other physi- 
cal contextsH3a. 

After introducing the order parameters in Section ||, 
we write down and classify quantum field theories for 
the critical points in Section III. We pause briefly in 



Section IV to consider the conventional transitions men- 



tioned above, in which the fermion bilinear order param- 
eter is not damped by Fermi surface excitations because 
of restrictions arising from momentum conservation, and 
show that it is described by field theories similar to those 
in Section [II. The latter are subjected to a renormaliza- 
tion group analysis in Section and the physical impli- 
cations of the results are noted in Section VI. 



II. ORDER PARAMETER AND PHYSICAL 
MOTIVATION 

Consider a correlated electronic system in two dimen- 
sions with enhanced SDW correlations at the wavevectors 
K x and K y ; these are vectors of equal length pointing 
along the x and y axes respectively. Similar considera- 
tions apply to other choices for the ordering wavevectors, 
but we will focus on this case because it directly applies 
to the doped cupratesES at doping concentrations above 
0.055. We can write for the spin operator S a (r, t) (r is 
a spatial co-ordinate, a — x,y, z labels spin components, 
and t is imaginary time): 

S„(r,r) = Re [e iK ^ xa (r, r) + e iK »- r $„ Q (r, r)] ; 

(2.1) 

Q xa and $> ya are then the SDW order parameters. Ex- 
cept of the case of two sublattice order (K^ = (n, 0) 
or (ir, 7r)), these order parameters are complex numbers. 
Concommitant with this SDW correlations, symmetry 
demandscJ'O that there are also "bond order" correla- 
tions at half the wavelength: 



Q a (r,T) ee S a (r,T)S a (r + a,r) 



Re 



.,21k. 



iK, 



(2.2) 



there is an implied summation over the repeated spin in- 
dex a. The vector a represents a bond (say the nearest 
neighbor vector), we have assumed that <& xa does not 
vary significantly over the spatial distance a, and the el- 
lipses denote numerous other similar terms which can be 
deduced from (2.1). For a a nearest-neighbor vector, Q a 



measures the modulations in the exchange energy, while 
for a = 0, Q a measures the local charge density (for the 
t-J model the charge density is linearly related to the on- 
site S^). So, quite generally, the bond order parameters, 
ya , measure modulations in all observables invariant 
under spin rotations and time reversal (such the electron 
kinetic, pairing, or exchange energies). The modulation 
in the total charge density may be strongly suppressed by 
the long-range Coulomb interactions, but this suppres- 
sion is not expected to apply to observables with a^O. 
We will loosely refer to such bond order modulations as 
CDW order in this paper, but the caveats discussed in 
this paragraph must be kept in mind. 

Now examine the structure of the effective potential, 
V(& xa ) controlling fluctuations of $ m (parallel consid- 
erations apply implicitly in the remainder of this..si.ibsec- 
tion to $j, Q ). On general symmetry groundsaE3, this 
potential has the form 



V(* xa ) = SKA 



-<b* cf>* <j> „<b 

2 *ia*»a* s 



x/3^x(3 



.(2.3) 



where s, u, v are phcnomcnological Landau parameters, 
and there is again an implied summation over repeated 
spin indices a,/3. In the usual Landau theory, the transi- 
tion to the onset long-range SDW order happens when s 
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first becomes negative, and the minimum of the potential 
moves from <fr xa = to a non-zero value. The theory of 
the critical point then accounts for fluctuations of $ m 
about the value of $ m = 0, and these allow for strong 
fluctuations in both the amplitudes and the phases of the 
three complex fields $ xa . 

In much of the discussion on stripe physics in the litera- 
ture (see e.g. Ref. ^ij), the physical picture advocated for 
the onset of long range SDW/CDW (or "stripe" ) order is 
quite different. As one decreases the doping from a Fermi 
liquid towards a Mott insulator, initial local stripe-like 
correlations form, the the orientation and phase of these 
regions changes from point to point; only when these lo- 
cally ordered regions align with each other, is there long 
range SDW/CDW order. In terms of our order param- 
eter formulation, this means that locally the amplitude 
$* a $ m first becomes large, and that long range order 
develops subsequently upon alignment of the orientation 
and phases of <& xa . 

In the familiar <p A field theory of a real vector order 
parameter, this change in perspective from a "soft-spin" 
Landau theory (with strong amplitude fluctuations) to 
a "hard-spin" perspective (with amplitude fluctuations 
quenched) does not make a fundamental difference in 
the physics being considered. The orientation fluctua- 
tions in the hard-spin approach are described by a non- 
linear a model, whose critical properties can be ana- 
lyzed by expansion away the lower critical dimension. 
Conversely, the soft-spin approach leads to an expansion 
away from the upper critical dimension. However, these 
are merely complementary approaches to the same criti- 
cal point, and the quantum numbers of the excitations in 
the phases flanking the critical point are the same in both 
approaches; indeed a l/N expansion (N is the number of 
field components) can interpolate between the soft-spin 
and hard-spin theories, aad also between the upper and 
lower critical dimensionsEa. 

For the present situation with a complex vector or- 
der parameter the situation is dramatically different, and 
there is a fundamental difference between the soft-spin 
and hard-spin app roac hes to fluctuations implied by the 
effective potential (2.3). The key difference is that not all 
order parameter configurations with a fixed &* a Q xa are 
physically equivalent, and they cannot all be rotated into 
each other by a symmetry transformation. A minimiza- 
tion of ( |2.3| ) at fixed $* a <& xa shows two distinct classes 
of minima, chosen by the sign of v; their most general 
form is: 



(I) v > : $ xa = n la + in 2a 

with 711,20 real, n\ a = n 2a and ni a n 2a = 
(II) zT < : $ xa = e i6x n xa 

with n xa real (2-4) 



By inserting into ( 2A ) , it is easy to see tha t (I) describes 
a spiral SDW for which the charge order in (2.2) vanishes, 



tant charge order. All experimental indicationsSB are 
that the spin fluctuations in the cuprates are always of 
the form (II) . It therefore seems physically reasonable to 
impose the amplitude constraint implied by (II) at an 
early stage, and to examine the theory associated with 
F($ la ) in the limit of large and negative v (stability re- 
quires u + v > 0). This last constraint is the more formal 
statement of the physical requirement that the spin and 
charge ordering is "stripe-like" . 

As has been discussed in Ref. |2^, the imposition of 
the limit v < at the outset leads to a theory for phase 
and orientation fluctuations with a non-trivial structure 
that is very naturally described by a Z 2 gauge theory. 
With the parameterization (II) in ( |2.4| ), and for fixed 
$* a $ la , the order parameter § xa takes values in the 
space (S 2 x S\)/Z 2l where the S 2 pertains to the orien- 
tation of n xa , the S\ is the phase factored e 19 * , while the 
all-important Z 2 quotient accounts for the invariance of 
<fr xa under n xa -^j— r-n xa and 9 X — ► 9 X + tt. We can de- 
fine global actionsEj for single-valued fields n xa and e lS:E 
only at the cost of introducing a Z 2 gauge field <7y = ±1. 



A simple phenore 
such an action ist 



■Jiological form for the spatial terms in 



Sz 2 = — - ^ cTjj cos(0 xi - 6 xj ) - — cr 

9s 1 ., 9n . . .. 

(w w) 



(2.5) 



□ □ 



while (II) is a collinear SDW with a modulation in the 
length of the spin as a function of r, and a concommi- 



where i, j are lattice sites, and gg, g n , and K are coupling 
constants that vary with doping. The Z 2 gauge invari- 
ance ensures that all physical properties are invariant un- 
der the Z 2 gauge transformation 9 xi — > 6 xi + (1 — o _ J )7r/2, 
n xa i — > n xa i<Ji and <xy — > o~i<JijO~j. The point defects as- 
sociated with order parameter space (£2 x S\)/Z 2 were 
analyzed in Ref. ^3|, and a fundamental role is played 
by the "stripe dislocation" , sketched in Fig [j], which is 
a half- vortex in the angular field X . This vortex also 
carries a Z 2 gauge flux of — 1 , as is easily seen from (^1]) . 

We are now in a position to introduce our candidate 
transition for strongly coupled quantum criticality in the 
presence of a Fermi surface. Begin in a Fermi liquid state, 
with only short range spin a nd charge order correlations: 
(®xa) = 0, ($l a ) = (s in ( |2~3l ) is large and positive, gg, 
g n in (2.5) are not large) and only a small core energy 
for the stripe dislocations, which therefore proliferate. 
The dislocation core energy is controlled by the Maxwell 
term for the Z 2 gauge field, K, and a small value of K 
implies that the Z 2 gauge field is in its confining state. In 
this state, the strong fluctuations of the Z 2 gauge field 
"bind" the n xa and e * fluctuations together, and the 
appropriate collective exci tatio n is the conventional field 
Q xa itself, as it appears in (2.1). The single $ m collective 
mo de therefore describes both the SDW fluctautions (via 
2T|)) and the bond order/CDW fluctuations (via (JT|)). 

Now increase the dislocation core energy by increasing 
K, while keeping J reasonably large (or gg, g n small): this 
will decrease the number of stripe dislocations but still 
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FIG. 1: A stripe disloc ation : gra y sc ale plot of 5p(r) ~ 
cos(2K ;c ■ r + 29 x ) (see (Q and (Q) with 0* = tt/2 + 
(1/2) tan _1 (y/:r) containing a half vortex. As noted below 
( j2.2j ) <5p(r) should not be literally interpreted as the electron 
charge density, and the modulation about could be associ- 
ated with e.g. the exchange per link. The quantum phase 
transitions discussed in this paper involve a transition from 
a Fermi liquid state in which these dislocations proliferate to 
one in which they are suppressed; there is no long-range spin 
or charge density wave order in either state (see also Fig. bj). 



maintain ($> xct ) = 0, ($j a ) = 0. In the cuprates we imag- 
ine that the increase in K is associated with the decrease 
in carrier concentration towards the Mott insulator. At 
a critical value of K c , the Z 2 gauge theory undergoes 
a transition into a deconfmed state in which the stripe 
dislocations are suppressed. The transition at K = K c 
is the focus of interest in this paper. In this deconfmed 
state, the n xa and e lBx fluctuations become independent 
collective excitations: i.e. the single spin/charge collec- 
tive mode, has frac tion alized into two independent 
modes. The bond order, ( |2.2| ), fluctuations are now given 
by the 'square' of the ejjf collective mode, while the spin 
fluctuations, (2.1), are the product of the n xa and e x 
collective modes. 

We reiterate that there is no long range spin or charge 
density wave order on either side of this fractionalization 
transition, and it is entirely associated with the suppres- 
sion of the defects in Fig M see Fig ||. At dopings lower 
than those shown in Fig H the cuprates possess phases 
with a variety of possible long range orders (e.g. SDW or 
CDW) which were reviewed in Ref. These transitions 
are associated with conventional order parameters, and 
we will not consider them in this paper, as theories are 
already available in the literature. 



III. QUANTUM FIELD THEORIES 

We are interested in describing the universal proper- 
ties of the transition at K = K c . These turn out to be 



<*a>=fJ 
<*a>=0 

< ty Xi y> = 

Deconfined phase of 
Z 2 gauge theory; 
Dislocations 
suppressed 



<*a>=0 

<<>=o 

Confining phase of 
Z 2 gauge theory; 
Dislocations 
proliferate 



K; 



K 



1 



FIG. 2: The quantum phase transition at K = K c is the one 
studied in this paper. The phases on both sides of K c have 
well defined electronic quasiparticles near a Fermi surface and 
no long range spin or charge order ((& x ,ya) = 0, ($ 1 ,,„ ) = 0). 



are defined in Section III: these 



The order parameters 
are Ising order parameters dual to the gauge theory. The 
phase on the right is a conventional Fermi liquid; that on the 
left is a Fermi liquid with local stripe correlations with ran- 
dom phases and spin orientations, and a 'topological rigidity' 
associated with the suppression of stripe dislocations. In the 
application to the cuprates, the horizontal axis represents in- 
creasing doping and the critical point is near optimal doping. 
States with long range spin or charge density wave order do 
appear at smaller doping but these have not been shown. The 
phase transition at K = K c should be distinguished from the 
conventional SDW ordering transition in a superconductor; 
thpsi|e«curs at lower doping and was used recently by one of 
usE3Ej to predict and explain the results of neutron scattering 
and STM experiments. 



controlled by two distinct characteristics of the electronic 
system. 

(?) Tetragonal or orthorhombic crystal symmetry: With 
tetragonal symmetry we have to simultaneously consider 
the separate Z 2 gauge theories associated with fraction- 
alization transitions in both $ m and & ya . With or- 
thorhombic symmetry, only one of these will be selected 
(the "stripes" have a preferred direction) . 
(m) Presence or absence of long-range Coulomb interac- 
tions: As we argue below, an important coupling of the 
critical Z 2 gauge degrees of freedom to the fermions arises 
from long-wavelength deformations of the shape of the 
Fermi surface. As is well known, in the presence of long- 
range Coulomb interactions the "dilatational" mode, in- 
volving local changes in the area of the Fermi surface, 
is strongly suppressed. For completeness we will also 
consider the case of fermions with only short-range in- 
teractions (perhaps the Coulomb interactions have been 
screened by a metallic gate), where the coupling to the 
dilatational mode has to be included. 

With these two criteria, four distinct categories 
emerge, which we label A, B, C, D, as summarized in 
Table Q. We will describe these cases in the subsections 
below. 
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TABLE I: Classification of quantum field theories. 



the bilinears in ( pj| ) and ( |2.2| ): 



Case 


Conditions 


Fields 


Action 


A 


Orthorhomic symmetry and 
long-range Coulomb interactions 






B 


Tetragonal symmetry and 
long-range Coulomb interactions 




Otii ~~t~ O T 
xy i l 


C 


Orthorhombic symmetry and 
short-range interactions 


4>x, i> 


Sx + s$ 


D 


Tetragonal symmetry and 
short-range interactions 


4>x, <t> v , 


Sxy + Si 



-AxX;cJ(r)^<v(r)S a (r) - A 2 £>t(r)<v (r)Q B (r), 

r r 

where r Q are the spin Pauli matrices. We now imagine 
integrating out $ M and generating terms which couple 
the electrons to the Z2 gauge degrees of freedom. The 
electrons are gauge singlets, and so can only couple to the 
physical Z2 gauge flux: the simplest such coupling is one 
in which the coupling K in ( |2.q ) becomes r, r dependent 
via its dependence on the local fermion bilinears 



A. Case A 

We begin with case A, for which the critical theory 
has the simplest form. As we noted in Section [n], the 
transition at K — K c is described by a Z2 gauge theory. 
We know that the pure 2+1 dimensionaL&i|gauge theory 
is dual a 2+1 dimensional Ising mode£jc3c3; we repre- 
sent this Ising order parameter by a real scalar field <j> x . 
For now, we neglect the independent Z2 gauge theory as- 
sociated with & ya , because orthorhombic symmetry has 
oriented all the spin and charge density waves along f la . 
The action for <p x in the vicinity of the transition is the 
familiar field theory for the Ising critical point: 



1 



d 2 rdr 



(3.1) 

we have chosen the spatial units so that the velocity of 
cj) x excitations is unity in both directions, u is a quartic 
non-linearity, and the quadratic coupling s ~ K — K c is 
the tuning parameter across the transition. The field <p x 
can be viewed as the creation/annihilation operator for 
the dislocations: (j> x is condensed in the confining phase 
where the dislocations proliferate, and (j> x has short-range 
correlations in the deconfined phase where dislocations 
are suppressed. 

In principle, we also need to account for the "matter" 
fields 9 X , n xa in (2.5) near this confinement transition. 
Closely related [ Za gauge models have been considered 
in the literatureL^ISa, and it has been established that 
their fluctuations do not change the universality class of 
the transition. Their presence does modify the effective 
couplings in ([O]), but the confinement transition con- 
tinues to be described by the Ising field theory (3.1) all 
along the phase boundary. 

The main question we need to address here is the cou- 
pling between <fr and the electronic excitations near the 
Fermi surface; the gapless excitations near the Fermi sur- 
face have the potential of inducing long-range interac- 
tions which are more destructive than the matter fields in 
The electrons do not carry any Ising gauge charges, 



{K, s) - (K, s) + A 3 ct(r) C? (r) + A 4 4(r)V^(r) + . . . , 

(3.3) 

where the ellipses represent terms with higher deriva- 
tives; we have to include all such terms because the typi- 
cal electron momentum is of order the Fermi momentum, 
and this is not small. Upon performing the duality tran- 
sition to the Ising field <p x , the mapping applies to s, the 
co-efficient of <f) x , via its dependence upon K, and this 
is also noted in (3.3). We can now integrate over the 



fermionic degrees of freedom and examine the structure 
of the terms which modify S x : it is not difficult to see 
that apart from the mod ification of the numerical values 
of the couplings in (|]l]), all such terms are formally ir- 
relevant. The leading non- analytic terms induced by the 
low energy fermionic excitations arise from long wave- 
length deformations of the Fermi surface: in the presence 
of long-range Coulomb interactions, all such excitations 
allowed by symmetry to couple to 4> x are suppressed. So 
we reach the important conclusion that the critical the- 
ory for case A is specified by (3.1). 

There is one important caveat to this conclusion that 
deserves highlighting. Our analysis has neglected topo- 
logical Berry phase terms that could be present in the 
coupling between (j> x , the fermions, and the Z2 gauge 
field (Tjj . Similar phase factors do arise in other analyses 
of Z 2 gauge theoriestS. We leave the study of this issue to 
future work, and restrict our attention here to couplings 
associated with (3.3). 



(2-5) 



and couple to degrees of freedom of interest here only via 



B. Case B 



The higher tetragonal symmetry has two important 
consequences: we have two consider fractionalization in 
both $ m and (fr ya simultaneously, and distortions of the 
Fermi surface with B\ symmetry do couple to the critical 
degrees of the freedom even in the presence of long-range 
Coulomb interactions. 

There are now two dual Ising fields, <f> x and <p y , asso- 
ciated with the independent Z2 gauge fields in the x and 
y directions; simple symmetry considerations show that 
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S x in (3.1) is replaced by 



S xy —Id rdr 



- {d T <t> x ) 2 + i (d T (j> y ) 2 + i {V r <j) x f 



1 



(3.4) 



for simplicity, we have neglected a possible difference 
in the velocity of the <f> x modes in the x and y direc- 
tions. The [Critical properties of the theory S xy are well 
understoodE.3: they are described by a stable fixed point 
with v* = and a global 0(2) (or XY) symmetry of 
rotations between <f> x . y . 

The fundamental new phenomenon here is that low 
energy deformations of the Fermi surface destabilize the 
critical point just noted. Although long-range Coulomb 
interactions do suppress a local uniform dilatational 
change in the Fermi surface, they do not suppress shape 
deformations which preserve the local volume. In partic- 
ular, the tetragonal symmetry allows an oscillation of the 
Fermi surface with B± (or d x i_ y -2) symmetry which has 
singular contributions at low energy, as has been empha- 
sized recently by Oganesyan et alrB; moreover, by sym- 
metry, this deformation will couple linearly to <p x — <p~. 
Integrating the-Germi surface deformations out by a stan- 
dard procedurealI3, we obtain the following non-analytic 
term in the effective action: 



f d 2 q 


f dto 


M | 


1 (2tt)2 J 


' 2^ 


|q| 1 



7T~ 777 hPxiq^J 



(3.5) 



the |w|/|q| factor is the standard density of states of 
long-wavelength Fermi surface oscillations in dimensions 
> 2. A simple scaling argument now gives the scaling di- 
mension of the coupling A5 at the v* = fix ed point 
of S xy . The |w|/|q| factor implies that (B^) couples 
(<^(r, r) — <f>y(r, r)) to itself at distinct spacetime points, 
which precludes any additional operator renormalization; 
consequently 



dim[A 5 ] =d + l- 2dim[4>l - <j> 2 } 



(3.6) 



here and henceforth we will state such scaling dimensions 
for general spatial dimensions d, although our interest is 
in d = 2. The scaling dimension of <j) x — (f> 2 c, 
be related to the anisotropy crossover exponent 4 



turn 

4>T 



dim[^-#] = d- 



(3.7) 



where v is the correlation length exponent of the v*. -p. ft 
fixed point of S xy . Using the known numerical valuescj'O 
of these exponents in d = 2, <f>T = 1.184 and v = 0.67155, 
we obtain dim[A 5 ] = 0.52 > 0. So we have established our 
earlier claim that B\ deformations of the Fermi surface 
constitute a relevant perturbation. 



It is now necessary to extend the renormalization group 
to include contributions higher order in A5. However, a 
simple oneJoop analysis shows that this is not a stable 
procedure! 4 !!: at order A§ we generate additional quar- 
tic term s wh ich have a different uj, q dependence than 
that in (3.5), and new relevant terms are generated at 
each successive order. The solution to this conundrum is 
to introduce a new real scalar field, tpi 1 as a Hubbard- 
Stratonovich decoupling of (3.5): this field will be a mea- 
sure of the amplitude of B\ deformations of the Fermi 
surface. Note that we expect A5 > 0, and so to make 
the Hubbard-Stratonovich transformation positive defi- 
nite, it is necessary to subtract a (4> x — (j) y ) 2 term before 
decoupling, and to compensate for this by corresponding 
changes to u,v. In this manner we get the action 



Sr 



d 2 q 

_ 9 
2 



du> 
2^ 



1 + o/f-r \ipi(q,u 



d'rdr^j {eg - <%) 



(3.8) 



our final effective action for case B is S xy +Sj. In obtain- 
ing Si we have expanded the coefficient of ip 2 for small 
A5 ~ 05, and rescaled ipi to make the co-efficient of the 
analytic ip 2 term unity. We now have two new coupling 
constants, g and etj, and it can be checked that A5 > 
implies a/ > 0. Further, at the fixed points of ^j,, we 
will show in Section VA that dim[aj] = dimfAs]. 

The most important property of S xy + Si is that this 
action is stable under renormalization, and no new rel- 
evant terms are generated at any order. Consequently, 
it can be subjected to a standard field-theoretic renor- 
malization group analysis, and this will be presented in 
Section 0. 



C. Case C 

The analysis of the remaining cases is a s traigh tforward 
generalization of the discussion in Section III B . 

As in Case A, the orthorhombic symmetry implies that 
we need only one Ising field <f) x . However, because of the 
absence of long-range Coulomb interactions, dilatational 
changes in the local Fermi surface volume with full lattice 
symmetry do carry large low energy spectral weight. We 
measure these changes by the real scalar field ip, and 
obtain its effective action 



Sit 



d 2 q 

_ 9 
2 



2^ 



1 + off h%,o;)r 



d 2 rdTip(j} 2 . 



The effective action for case C is now iS : 
gument similar to that in Section 



IIIB 



(3.9) 

S,p. An ar- 
can be used to 



compute the scaling dimension of a: as in (3.6,3.7) we 
find that the fixed point describing the transition of S x 



dim [a] 



1 



a 
v 



(3.10) 
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where a, v are the standard exponents of the 24=1 di- 
mensional Ising model. Using their known valuescj, we 
obtain dim[a] = 0.174 > 0. So the fixed point of case A is 
unstable to a Fermi surface with short-range interactions, 
and there is new universality class of critical behavior. 

We note that the discussion in this subsection, with 
zero sound modes of a Fermi surface coupling to the 
square of an order parameter, bears a passing resem- 
blance to a discussion of the influence of phonon modes 
in the superfluid state of bosons presented by Frey and 
BalentsH 



as the composite fields (j) x y did in Section III . For or- 
thorhombic symmetry and short-range interactions (as 
in Case C), only $ m is critical (say), and its coupling to 
the long-wavelength particle-hole continuum leads to the 
term 



A 6 



<Pq 
(2nf 



dlM \uj\ | 2 ,2 



which is the analog of (3.5). Just as in 



dim[A 6 ] = 



(4.1) 



(4.2) 



D. Case D 

With tetragonal symmetry and no long-range Coulomb 
interactions, both dilatational and B\ modifications of 
the Fermi surface have to be included. So the full effective 



action is S xy + Si 



SL where 



SL — — 



d 2 q 



5" 



— / d 2 rdrip 



At the 



fixed point of S xy , we can show, as in 



Section III C , that dim [a] = a/v, where a, v are now 
critical exponents of the 2+1 dimensional XY model. 
Numerically^ we have dim [a] = —0.022 < 0, and so a is 
irrelevant at this fixed point. However, the perturbations 
in Si remain relevant, and we cannot a priori conclude 
that the terms in SL will also be irrelevant at any pos- 
sible sta ble fi xed point of S xy + Si; so it is necessary to 
include (3.11) in the analysis of Case D. 



and so the coupling in ( |4.l[ ) is relevant if and only if 
a > 0. Most conventional order parameters have a < 
in 2+1 dimensions, in which case the Fermi surface has no 
influence on the critical theory. The main exception is the 
Ising case, found e.g. at a two sublattice CDW ordering 
transition, in wh ich case the theory will identical to that 



in Section III C 



For a system with tetragonal symmetry, there is 
stronger coupling to quadrupolar distortions of the Fermi 
surfa ce, just as discussed in Section IIIB. The analog of 
(3.11) (|J) is the term 



1 



{2-Kf 



dw \lo\ 
2^N 



(q,w) - \§ ya \ (q,w) 



(4.3) 

and the scaling dimension of A7 is related to an anisotropy 
crossover exponent; this always such that A7 is strongly 
relevant. The analysis of the resulting theory will be 
similar to Case B above, although we will not present it 
explicitly. 



IV. CONVENTIONAL ORDER PARAMETERS 
WITHOUT FERMI SURFACE DAMPING 

This section makes a brief detour to a separate class of 
quantum phase transitions mentioned in Section |: those 
with conventional order parameters which do not couple 
linearly to low energy excitations on the Fermi surface. 
We will consider here a simple example of such a tran- 
sition, although the results discussed here are far more 
general: an SDW orderi ng t ransition with the order pa- 
rameters defined as in (2.1), from a Fermi liquid state 
with (<& x .ti« 



to an SDW state with ($ 2 



^ 0. 



Further, we assume that the values of K. x ,y are such that 
they cannot connect any two points on the Fermi surface; 
consequently the <fr x ^ ya fluctuations are undamped and 
the quantum critical modes are strongly coupled. The 
effective action for Q x , ya will have a structure similar to 



that of S x , S xy in (3.1), (3.4), apart from ch ang es in the 
quartic terms which are similar to those in (2.3). 

It should now be clear that the zero momentum com- 
posite fields l^z.yal 2 can couple to long-wavelength de- 
formations of the Fermi surface in much the same way 



V. RENORMALIZATION GROUP ANALYSIS 

This section returns to the collective mode fraction- 
alization transitions considered in Sections |n] and III; 
closely rel ated analyses will apply to the cases mentioned 
in Section IV 

Case A is described by the purely bosonic theory ( |3.1| ) , 
whose properties are very well known: its critical point 
is in the universality class of the 2+1 dimensional Ising 
model. So we move directly to case B in the first subsec- 
tion, and mention the straightforward generalization to 
the other cases later. 



Case B 



Renormalization group (RG) analyses of models with 
a structure similar to S xy + Si have appeared elsewhere 
in the literature (see e.g. Ref. [l9|), and so we can omit 
most of the technical details. We consider the rescaling 
transformation 



r e 



t e 



(5.1) 
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under which the fields 4> x ,y, ipi transform as 



Mr,r) = e -( d+z+ W 2 V^(r',r'); 



(5.2) 



here z, 6, and rj are ^-dependent parameters to be deter- 
mined by loop corrections to the effective action under 
the RG transformation — the values of these parameters 
at any fixed point of the RG equations will specify the 
corresponding critical exponents. The one- loop graphs 
have the same topology as those considered in Ref. 
and we defer details on the evaluation of the loop inte- 
grals to Appendix |b|; the results are 
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(l-z-6)ai 

{A-d-z-9- 2rj) (u v\ 

2 9 ~ { 3 _ 2 J 9l{ °> Ql) 

g 3 l(l, ai ) 

(I0u 2 + 6uv + 9v 2 ) ^ n 
(4 — a — z — 2r\)u 1(0, a/) 

6ug 2 l(l,a T ) - 6.g 4 X(2, a/ ) 

(4-d-z-2r 1 )v- ( 4UV+ 2 3v2) l(0,a I ) 

2(4u + 3v) 2 

g 2(1, a/) 



-5 2 ^(0,a/) 

5 2 ^(«/) 

g 2 JC(a!) - r] 



(5.3) 



where I(m, a/) (m integer), J7(a/), and/C(a/) are scheme 
dependent functions which are specified in Appendix 
We have assumed that the strongly relevant "mass" term 
s has been tuned to its critical value, and are examining 
above only the flow with in t his critical manifold. The 
first of the equations in (5.3) is the most important: it 
is an exact equation (although the values of z and 8 do 
have corrections at higher order) which follows from the 
fact that there are no loop corrections to the non-analytic 
term with coefficient aj in the effective action. 

First, let us examine ( |5.3] ) in the subspace with aj = 0, 
where the long-range forces associated with the Fermi 
surface are absent. The results in Appendix 5.3) show 
that T(m, 0) = C independent of m, and J(0) = /C(0) = 
0; here C is a scheme-dependent constant whose value 
does not effect the physica l pr operties. It is then easy to 
search for fixed points of ( |5.3| ): there are 8 fixed points, 
and all of them obey 2 = 1 (this is exact to all orders for 



where e = 3 — d. The physical meaning of these fixed 
points becomes clear after integrating out the -0/ field in 
S xy at a j = 0: for this value of aj this yields an effective 
action with the same structure as S xy in ( |3.4[ ) but with 
modified coupling constants 



u 
v 



u-3g 
v + 2g 2 



(5.5) 



The reader can now observe that the fixed points in (5.4) 
all separate into pairs with equal values of u and v, and 
our notation for the fixed points has been chosen to make 
this evident: G\ and G 2 form one such pair, and similarly 
for the remainder. 

At this staffi©i we can compare results with the analysis 
of Aharonyrli 50 ! of the consequences of a cubic anisotropy 
on fixed points with 0(n) symmetry: he analyzed the 
analog of the model S xy for n-component fields, but with 
a somewhat different notation. The fixed points (5.4) 
are in complete agreement with his results: Gi,2 are the 
Gaussian fixed points, are the "cubic" fixed points, 
XY12 are the XY fixed points with 0(2) symmetry and 
no cubic anisotropy, and 1% % are the Ising fixed points 
at which the fields 0i. 2 are decoupled. As shown by 
Aharony, an XY fixed point is stable for the case of a 
two-component field being considered here. 

Here, we have to perform some further analysis to dis- 
tinguish between the two possible XY fixed points: these 
two have the same values of u, v, but only one has g* ^= 0. 
A standard computation shows that fixed point XY2 with 
g* ^ is the only stable fixed point among all those in 
( |5.4| ). The theory prefers the fixed point with g* ^= be- 
cause the fluctuations of ipi ~ cf> 2 — (j)?, are characterized 
by a divergent susceptibility at the critical pointEj. 

We are now in a position to consider non-zero values 
of a/. For small ai, the above analysis shows that near 
the fixed point XY? 



dai 
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3e 



-ar, 



(5.6) 
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FIG. 3: Plot of (1 - z- 8)/g 2 as determined by (Q. There 
is no zero crossing, and hence no finite, positive, non-zero ai 
fixed point. The full line is the minimal subtraction scheme, 
while the dashed line is the fixed dimension scheme. 



In other words, ai is a relevant perturbation with 
dim[a/] = 3e/5. We argued in Section III B that 
dim[a/] = 24>t/v — d— 1, and it can be checked that this 
is consistent with our present result and the e expansion 
resulted for 4>t/v. 

For general a/ , we have to consider the complete flow 
equations in (5J3). While these are quite complicated, it 
is clear from the first flow equation that a finite and non- 
zero ai fixed point, ai = a*j can only exist if and only if 
there is a solution to 1 — z — 8 = 0. From (|j^) we observe 
that at the one loop order, this condition translates into 
an equation for a J alone: 



FIG. 4: An in Fig H, but for (5.£) 



21(0,a* I ) + J(a})-IC(a* I 



0. 



(5.7) 



Using the results in Appendix H we plot half the left 
hand side of this equation in Fig |^ for two different 
renormalization schemes. No solution is found in either 
case, or in the additional cutoff schemes discussed in Ap- 
pendix |^. The value of (1 — z — 9) does initially decrease 
from its value at ai — which led to ([T(]), but this de- 
crease is not strong enough to lead to a zero crossing. So 
we have not found a finite a/ fixed point, but we can- 
not rule out the possibility that one will appear upon 
including terms higher order in g, m, v. 

It is clear now that any non-zero, positive ai runs away 
to a j = oo, and we have to consider the theory in this 
limit. An examination of the limiting values of the re- 
sults in Appendix ^ showed that there were no finite 
fixed point values for u, v in the limit of large ai either. 
In principle, it was possible that both a/ and g became 
large in a manner that permitted u and v to approach 
a finite fixed point value. However simple analytic and 
numerical analyses of (5.3) showed that this was not the 
case, and we always had flows to strong coupling for u 
and v also. Starting close to the stable fixed point at 
a = 0, we found that the flows sent v — > oo and u — > — oo. 
This is possibly suggestive of a first-order transition, but 



no definitive conclusion can be reached because the RG 
equations themselves breakdown once u or v become of 
order unity. 



B. Case C 

The analysis and results for Case C are similar to those 
for Case B above, and so we w ill be brief. 
The RG flow equations ( |5.3|) are replaced by 



da 
~dt 
(lg_ 

df 

du 

d£ 



(l-z-e)a 
(4 - d - z- 



277) 



yl(0, a) 



3?/ 2 

2r ) )u-—l(0,a) 



,g 3 X(l,a) 

(4 - d - z 
6ug 2 l(l,a) - 6g 4 l(2,a) 



= 1 



g 2 JC{a) - V 



(5.8) 



where the functions X, J \ K, remain the ones specified 
in Appendix These equations do have fixed points at 
a — M, which have been specified completely in earlier 
workEi However, the stable fixed point for a — is 
unstable to a non-zero a, and the condition for a finite 
and non-zero a fixed point is (replacing (5.7)): 



1(0, a*) + J(a*) -tC(a*) = 0. 



(5.9) 



The left hand side of this equation is plotted in Fig. [|, 
and again there is no solution. We find instead flows to 
strong coupling, as in Case B. 

Finally, we note that the flow equations and results for 
Case D are a combination of those presented already for 
Cases C and D, and we will not present them explicitly. 
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VI. PHYSICAL PROPERTIES 

This section considers some simple observable physical 
properties of the fractionalization transition described in 
We will find that the properties of many 



Section III 



observables are quite similar to another case with an or- 
der parameter which is a composite of the electron spin 
operator: the onset of spin glass order. To help make 
this comparison, we review the physical properties of the 
quantum spin glass transition in Appendix [A|. 

Our analysis of collective mode fractionalization theo- 
ries has found only one stable fixed point: that for case 
A with orthorhombic symmetry and long-range Coulomb 
interactions. The other cases exhibited flows to strong 
coupling. This section will study the general scaling 
properties of various physical observables for case A, and 
also extrapolate the results to possible strongly coupling 
fixed points for the other cases. 



A. Dynamic spin susceptibility 

We focus on the spin susceptibility near the wavevec- 
tors K K j,; other susceptibilities of conventional order pa- 
rameters can be treated in a similar way. The order pa- 
rameters associated with these susceptibilities, $> XtVa , are 
controlled by an effective action of the form (w n is a Mat- 
subara frequency at a temperature T) 



T 



E 



9 



(T\u n \ + A<? + B) |* 



along with a similar term with (x — ► y). Unlike the situa- 
tion in Section IV , we are considering the case where K Ii!( 
do connect two points on the Fermi surface, and hence 
there is over-damping of <& x ,ya fluctuations represented 
by the term proportional to T. Also, the <& x ,yu fluctu- 
ations remain non-critical with (& x . ya ) — at the frac- 
tionalization transition, and hence the parameter B > 0. 

We now have to consider the coupling of & x . ya to the 
critical modes 4> x , y - As noted at the end of Section fll A|, 



we ignore possible non-local topological terms, and focus 
on the simplest local couplings allowed by symmetry; for 
the cases with tetragonal symmetry (cases B and D), 
these are 



d 2 rdr 



Ci(l*-| 2 + |$ att | 2 ) {4>l + 4> 2 y ) 

C 2 (|<f,a| 2 -|<IV| 2 ) (^"#), (6.2) 



while for orthorhombic symmetry we merely drop the 
non-critical 4> y field above. We will consider the influ- 
ence of the critical <j> x ,y modes on <S> x ^ ya in a perturbation 
theory in (j. For case A, where a stable fixed point was 
found, such a procedure appears reasonably safe. The va- 
lidity for the other cases, with flows to strong coupling, 
has a greater uncertainty and this should be kept in mind 
below. 



The influence of the critical fluctuations on the non- 
critical & Xtya susceptibility will be determined by two 
(and higher) point correlators of the operators <j% — cf) 2 
and (fy%.+4>y f° r tetragonal symmetry, and 4> 2 alone for or- 
thorhombic symmetry. Let us represent these operators 
generically by O and its scaling dimension by 



dim[C] = do- 



(6.3) 



For case A, we have only O = 4> 2 and do = d + 1 — 1/f, 
where v is the exponent of 2+1 dimensional Ising model; 
for other cases, the value of do will be determined by 
the hypothetical strong-coupling fixed point. In pertur- 
bation theory in (^2 we can write for the dynamic spin 
susceptibility 



1 



-iTuj + Aq 2 + B + M (q, u) 



(6.4) 



with the "self energy" 



-q, e n +uj n ) 



M(q,w„) = C 2 T^ J ^^D (p,e n )x*(p 

(6.5) 

where £ is related to £1,2, and Do is the two-point cor- 
relator of O. We are intereste d in the behavior of M at 
small q and w; after inserting (6.4) into (HI]), we use the 



fact that B > to expand the momentum and frequency 
dependence of x$ in powers of 1/B. For the q depen- 
dence of M, a simple power-counting analysis shows that 
all terms in this series are free of infrared divergences 
as long as do > 0, in which case the q dependence of 
M is analytic and free of critical singularities. Case A 
easily satisfies do > 0, and we expect this is also true 
for the other cases, because a negative do would imply 
that correlations of O increase with increasing spacetime 
distance. 

In contrast to its q dependence, the to dependence of 
M is singular even for do > 0. This is most easily com- 
puted by analytically continuing (6.5) to real frequencies 
and taking its imaginary part. A standard analysis then 
shows that M(q,u>) has a q independent singular contri- 
bution in the quantum critical region which obeys the 
scaling form 



M(q,w) ~ T( 2do+z )/*Fi(w/T) 



(6.6) 



(we are using H = ks = 1), where F\ is a universal 
scaling function determined by the non-trivial interacting 
fixed point of the strongly- coupled quantum field theories 
considered in Section III; we expect Im.Fi(y — * 0) ~ y 
and Fi(y -> 00) ~ sgn(y)|?/|( 2cio+1) / z . Specifically, F x is 
related to the local correlator of 0, which has a singular 
contribution obeying the scaling form 



d 2 p 

(27T) 2 



D (p, LJn) - T 2d °/ z F 2 (iu n /T), (6.7) 

(the scaling function F 2 is a property of the strongly- 
interacting fixed point of the quantum field theories of 
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Section III) by the integral relation 



1m[Fi(y)} 

y dn 

— mm [F 2 (y - ft)} (1 + n(O) + n(y - ft)) 

n 7T 

— Qlm [F 2 (n + y)} (n(fi) - n(0 + y)) 

,) IT 

— 01m [F 2 (Q - y)\ - y) - n(fi)) , (6.8) 

- 1) is the Bose function at unit tem- 

arises 



where n(fl) — l/(e r 

perature. Note that the leading factor of in (6 



from the linearly de nsit y of states of low frequency exci- 
tations implied by (6.4), upon ig norin g the influence of 
M in \<s> on the right hand side of (6.5). This is permissi- 
ble as lon g as the contribution of M here is subdominant: 
from ( |6.6[ ) we see that this is so as long as do > 0. If this 
inequality is violated, then the perturbative expansion in 
£ fails, and a self consistent approach appea rs neces sary. 

It is interesting to note that the results (|3~i,6.6) im- 
ply a structure for the dynamic spin susceptibility which 
is closely related to phenomenological forjas that have 
been proposed recently by Schroder et al.cS to describe 
neutron scattering measurements on CeCug-^Aua;. Ele- 
gant theoretical models of a "local quantum criticality" 
based on the self-consistent treatment of quantum.impu- 
rity ( ' Ro ndo ' ) japdels have also been advancedEjOO and 
discussedcZlGaO'LZl as explanations for the observations of 
Schroder et at. Some-important differences between the 
theories of Si et alBO and our work here should be 
noted: underlying our results is a bona fide d + 1 dimen- 



sional quantum field theory (presented in Section III), 
and not a quantum impurity model. Our theory leads 
to a q-independent singular contribution to the dynamic 
spin susceptibility because B > 0, as shown above. The 
condition B > also implies the static spin susceptibility 
is not divergent, in contrast to impurity model analyses 
in which it does diverge at the impurity criticalpointEjE-X 
An associated fact is that the exponent in (S.6) is likely 
to satisfy (2do + z) / z > 1, although this need not be true 
for the spin glass transitions discussed in Appendix |Aj It 
has also been notecloEJO that the divergence of the lo- 
cal susceptibility could preempt .a_local treatment of the 
impurity model critical pointscJO. 

It should be clear that the above mechanism for 'local' 
quantum criticality in the dynamic spin susceptibility is 
quite a general property of a situation in which the pri- 
mary order parameter, which exhibits strongly coupled 
bulk criticality, is coupled to t he s pin operator only with 
non-linear terms like those in (6.2). As an illustration we 
recall in Appendix ^ that related considerations apply 
to the quantum spin glass transition. 



the electrons via couplings like 

C 3 / d 2 rdrO(r,r)4(r,T) C? (r,r) + 



(6.9) 



the ellipses represent additional terms with gradients of 
the electron operators (these are not small because the 
typical electron momentum resides on the Fermi surface). 

We should note at the outset that the status of a per- 
turbative computation of the electron self energy, S, re- 
mains open. For the case of overdamped fermion critical 
points below their upper critical dimension, the_corre- 
sponding issue was addressed by Altshuler et alE5 and 
Castellani et al!Ei who concluded that the leading per- 
turbative correction also gave the ultimate critical sin- 
gularity. The extension of this result to the strongly 
coupled critical points being discussed here is less clear. 
Indeed, for the case of a strongly coupled spin glass quan- 
tum critical point, the perturbative computation of the 
electron self energy is known to break downO (see also 
Appendix |A|) . 

With this caution in mind, we proceed with a com- 
putation of the consequences of (3.E). The electron self 



energy, S to order Q is 



1 



i(e n + uj n ) - e(p + q) - £(p + q, e„ + w n 



(6.10) 



where s(jp) is the electron dispersion. This expression can 
be evaluated using the methods discussed in e.g. Ref. |5^. 
The result depends upon the scaling function associated 
with Do for which we write 



£>o(q,w)~T( 



2d L 



-^/'^(w/T.q/T 1 /*), (6.11) 



with i<3 a non-trivial universal scaling function deter- 
mined by the strongly interacting fixed points of See- 
As in Ref. 



tion 



III 



the singular part of the self 
energy is momentum independent at this order, as long 
as z > 1. Its frequency and temperature dependence is 
specified by the scaling form 

E(q, w) - T^-W'Fifa/T), (6.12) 

with the scaling function F^ given by 

f°° dQ f 00 

Im[F 4 (y)] = — k d - 2 dk!m[F 3 (n,k)] 

JQ 7T Jo 

x (2n(fi) + f(n -y) + /(fi + y)) , (6.13) 



where n(fi) was defined below ( |6.8[ ) and f(Q) = l/(e n 
1) is the Fermi function at unit temperature. 



B. Fermion self energy 



C. Transport 



We follow the same perturbative approach discussed in 
Section VIA. Singular fluctuations of O will couple to 



We identify two distinct contributions to the critical 
electrical transport properties: the electronic quasiparti- 
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cle and the bosonic collective modes. However, the cau- 
tions about a perturbative treatment of the coupling be- 
tween these modes mentioned in Section VI B apply also 
to our discussion here. 

The electronic contribution is characterized by the 
transport relaxation time r tr . The electron scattering 



discussed in Section VI B all occurs at small momenta p 
and, as is well known, the transport rate then acquires an 
additional factor of p 2 over the quasiparticle scattering 
rate in the integral over scattere d mo menta. Extending 
the scaling arguments of Section VI B we see easily that 



1 



rp(2do + l)/z 



(6.14) 



The collective mode contribution is the analog of the 
"Aslamazov-Larkin" term in the theory of fluctuating 
superconductivity: fluctuations of the order parameters 
4> x ,y may themselves contribute to the electrical current 
J. The relationship between the dual fields <f> x ,y and the 
microscopic degrees of freedom is extremely complex, and 
so we will restrict ourselves here to deducing the expres- 
sions for the current operator using symmetry arguments 
alone. The current J is odd under time reversal and spa- 
tial inversion, while the fields <j) x . y are even under both. 
Furthermore, the current should be invariant under the 
dual Ising symmetries involving the change in sign of <f> x 
and/or <p y . This indicates that the current operator has 
contributions from operators like 



Jx ^ <9r 4*x &x 4>x , 4>xd T d x 4>x , d T 4>yd x 4>y 



(6.15) 



and similarly for J y . The scaling dimensions of these 
operators can be determined by standard methods for 
case A (some combination is related to components of 
the stress-energy tensor with scaling dimension d + 1, 
but the others require a (3 — d) expansion), while for 
the other cases the scaling dimensions remain unknown 
quantities to be determined by the hypothetical strong- 
coupling fixed point. If we gcncrically denote the scaling 
dimension of these operators by dj, the Kubo formula 
implies that they will yield a singular contribution to the 
conductivity, a, which scales as 



8a ~ T(?dj-d-2z)/z 



(6.16) 



Note that because the current is carried by both the elec- 
tronic and bosonic degrees of freedom, we expect that we 
do not have to worry about constraints imposed by over- 
all current conservation, as is the case in theories-of the 
superfluid-insulator transition in bosonic systems^. 



liquid state in which the dislocation in Fig [j] prolifer- 
ates, to a 'topologically ordered' Fermi liquid state in 
which these dislocations have been suppressed. Note 
that there is no long-range SDW/CDW order in cither 
state, although the amplitude of these orders must be lo- 
cally large to allow identification of the dislocations. The 
conventional SDW/CDW orders are suppressed in both 
phases by 'spin- wave' fluctuations and by the prolifera- 
tion of other defects, such as the hedgehogs in n xa and 
integer vortices in 9 X . We discussed some simple observ- 
able properties of such a topological phase transition in 
a Fermi liquid, and found results which were quite simi- 
lar to those near quantum spin glass transitions, and to 
the phenomenaLogical 'local' scaling forms proposed by 
Schroder et alSB. Despite their 'local' nature, the crit- 
ical singularities are controlled by a bulk quantum field 
theory, and not a quantum impurity model. 

The dislocation proliferation transition discussed here 
has been proposed as a candidate for a quantum 
critical point y p jj lr l optimal doping in the cuprate 
superconductorsacfl. In this context, it is interesting 



note that recent scanning tunnelling microscopy image: 
of the surface of Bi2Sr2CaCu20s+5 clearly show disloca- 
tions with the geometry of Fig [l] (see Fig 2c of Ref. 62). 

The main shortcoming of our analysis here has been 
runaway flow of the renormalization group flow to a non- 
perturbative regime for the interesting cases. It is pos- 
sible that this signals a more fundamental change in the 
nature of excitations near the Fermi surface than is ac- 
counted for by our formalism. We integrate out the low 
energy fermions at an early stage, and the strong cou- 
pling situation may benefit from an analysis in which 
quasiparticle degrees of freedom are retained more ex- 
plicitly. Future work in this direction, along the lines of 
Refs. |63L appears desirable. 
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APPENDIX A: QUANTUM SPIN GLASS 
TRANSITION 



VII. CONCLUSIONS 



This paper has described strongly interacting quan- 
tum phase transitions in the presence of a Fermi sur- 
face. The primary examples considered involve the phase 
transition described in Fig |^, from a conventional Fermi 



The quantum spin glass transition in metallic 
systeipSjJp^beeaj-studicd in a number of recent 
worksaoi£3'^oE30E3. Here we review some physical 
properties wh ose scaling structure resembles the results 



in Section VI 



The central actor in the theory of the spin glass tran- 
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sition is the order parameter functional 

Q^(r,ri,r a )~5S(r s T 1 )S|(r > 7i 1 ) (Al) 

where a, b are replica indices. In the quantum critical 
region, this functional has an expectation value of the 
form 

(Q^(r,n,r 2 )) = S a0 S ab D(n - r 2 ) (A2) 

If the quantum critical point obeys hyperscaling prop- 
erties, then Fourier transform of D(t) will obey "cj/T 
scaling" : 



D( W ) = W 5 ( W /T) 



(A3) 



It should be noted that the mean-field solution of the 
metallic spin glass does not obey (|A3|), and instead has 
a characteristic frequency scale ~ T 6/2 . However, other 
mean-field spin glass models which account for proximity 
to a Mott insulator do obey lu/T scalingBBBH; These 
latter solutions have the value /j, = 0, and D(lu) has the 
form 







" ln (r) 









2fi/(0) 
^ 2 - {u/TY 



ITT ( Ul \ 

— tanh — — , 
2 \2TJ 



(A4) 



where /(f2) was defined below (3.13) — for this value of 
fi there is an additive real contribution dependent upon 
the cutoff A, but the remainder is a universal function 
of lu/T, and the function D{lu) is analytic as lu — > 
at any nonzero T. We will assume below that the low- 
dimensional spin-glass quantum critical point does obey 
the hyperscaling property (A3), even if the mean-field 
theory for its particular situation does not. 

We now consider the computations of the physical 
quantities discussed in Section |vi[ For the dynamic spin 
susceptibility near the wavevector K, note that symme- 
try allows a coupling like 

C 4 J d 2 rdr 1 dr 2 g^(r,T 1 ,r 2 )Re [*£(r, nj^r, r 2 )] , 

(A5) 

along with the corresponding (x —> y). Using the ex- 
pectation value (A2) in (A5) , wc imme diat ely obtain a 
dynamic spin susceptibility of the form (fTJI) with 



M(q,w) ~ D(w). 



(A6) 



The singular part of the inverse susceptibility is q inde- 
pendent and a function of lu/T, as in the local quantum 
criticality scenario. Again, however, the scaling function 
in (A3) is determined by a bulk d-dimcnsional quantum 
field theory. 

Finally, we turn to the electron self energy, and the 
associated single particle and transport lifetimes. One 
approach, which follows Ref. [36[ is to compute this per- 
turbatively from the coupling which follows directly from 



the identification ( |Al| ) 

xc^(r,T 2 )r^,c b v ,(r,T 2 ). (A7) 



Using (A2), we can easily compute the electron self en- 
ergy, E perturbatively in £5, and at order £| we find the 
contribution 



E(q,w„) ~ ($Tj2sgn{e n )D(u; n + e n ) 



(A8) 



which is momentum independent, but a sing ular function 
of frequency and temperature. Using (A3), we deduce 
that E also obeys the scaling form 



E(q,w„)~ T^i^w/T) 



(A9) 



with 



Im[F 6 (y)] 



Im[F 5 (fi)] 



x (2n(fi) + f(Q -y) + /(fi + y)) , (A10) 



a relationship similar to ( 6.1 3| ) . Note again the local na- 
ture of the self energy. As the expression (A8) scatters 
the electron by large momenta around the Fermi sur- 
face, there is little difference between single-particle and 
transport scat terin g rates, and the latter can be deduced 
directly from (A9). It should also be noted here that the 
present perturbative computation for E is known to break 
down in certain cases: this was pointed out by Parcollct 
and Georged£2l who showed that the electron self energy 
was far more singular than ( |A9| ) for the spin glass models 
of Ref. 61 64|,|67] which had fi = 0. Nevertheless, the local 
nature of the self energy, and u>/T scaling was preserved. 



APPENDIX B: COMPUTATIONS FOR RG 
EQUATIONS 

This appendix will present some further details of the 
derivation of the renormalization group equations for 
Case B, described by the action S xy + Si. We noted 



in Section V A that the topological structure of the dia- 
grams is identical to those considered in Ref ^9|, and so 
we will not reproduce the diagrams here. 

The exponents r\ and z arise from the momentum and 
frequency dependence of the self energy E^ of the c^i. 2 
field. To leading order, the only such contribution to E 
is 



E (g,w) = g 



o / d d p dfl |p| 



(2ir) d 2vr (|p| + a/ |n|) 
1 

X ((p + q) 2 + (w + fi) 2 + s) 



(Bl) 
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where g 
we find 



|q|. By an expansion of for small g and to, where we have dropped the subscript Zona for brevity. 



X2 ) 



V + 2(z-l) = (-_)E M (0, W ) 



d 
dq 
d 

duj 



9=0 



(B2) 



j=0 



Here the s ubs cript A is meant to indicate that the in- 
tegrals in (Bl) are to be carried out as part of a renor- 
malization group procedure. The precise nature of the 
integral then depends upon the scheme being used: 
(i) in the dimensional regularization/minimal subtrac- 
tion scheme, we perform the integrals in arbitrary di- 
mension d, and pick out online the poles in e = 3 — d in 



(B2) 



(ii) in the fixed dimension approachO, we perform the 
integral in d = 2 for a finite s, take the derivative of the 
result with respect to *Js, and set_s.= 1. 
(Hi) in the soft-cutoff approachuL 3 ], we perform the in- 
tegrals with a smooth cutoff on the scale A, take the A 
derivative, and then set A = 1. 

Each of these approaches yields closely related expres- 
sions for J(ai) and K.(ai). 

The remaining terms in (|]^) arise from simple one- 
loop integrals which can be evaluated at zero external 
momenta and frequency; these yield the expression 



l(m, aj) 



d d p dn 



(2ir) d 2tt (|p| + ai\n\) m (p 2 + n 2 + s) 2 ' 

(B3) 

Again the A subscript indicates that the integrals are to 
be evaluated in o ne o f the renormalization group schemes 
discussed below (p2j). 

We can now present our explicit results for l(n,ai), 
J(a{), and /C(a/), for the respective schemes. 

In the dimensional regularization approach, picking 
poles in e = 3 — d, we obtained 



1(0, a) 
1(1, a) 

1(2, a) 



1 

8^ 

4a 3 In a 



2a 3 



3na 2 - 2a 



87r 3 (l + a 2 ) 2 
4a 5 - 3?ra 4 + 16a 3 In a + 6na 2 -4a + n 



(3 



J(a) 
K(a) = -3J(a) 



8tt 3 (1 + 
2 In a)a 5 — 37ra 4 



,2\3 



(2 + 6 In a)a 3 



127r 3 (l + a 2 ) 3 



(B4) 



Similar expressions were obtained in the fixed- 
dimension approach in d = 2: 



1(0, a) 
1(1, a) 

1(2, a) 

J(a) 
JC(a) 



1 

8^ 



(l-a)v / l + a 2 +a 2 £(q) 
87r(l + a 2 ) 3 / 2 



(a 3 - 2a 2 - 2a + 1) VTT^ 2 + 3a 2 £(a) 

87r(l + a 2 ) 5 / 2 
a(2a 2 - 3a - l)\/l + a 2 - a 2 (a 2 - 2)C(a) 



16tt(1 + a 2 ) 5 / 2 



-2.7(a) 



(B5) 



where 



C(a) = In 



(B6) 



We also obtained results in scheme (Hi) above, with 
both soft and hard cutoffs; however we refrain from giving 
complete expressions as the results are similar to those 
above. 
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